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Abstract 

The A decay in the nuclear medium is calculated in the relativistic meson- 
nucleon model. The delta spreading width is calculated and compared with 
the Pauli-blocked 7rN decay width. The influence of relativistic mean fields is 
also studied. We stress the importance of understanding the delta spreading 
width in interpreting experiments involving delta resonances. 

I. INTRODUCTION 

The understanding of the properties of the A resonance in the nuclear medium, especially 
the branching ratios of its various decay modes, is important for the analysis of scattering 
processes involving A-hole excitations. 

A A, which mainly decays into 7rN in free space, acquires additional decay channels in 
the medium. For example, it can decay through the mechanism, A + N — > N + N, commonly 
known as "delta spreading width" which can be seen as reabsorption of the virtual pion 
through a particle-hole excitation of the nucleus. Indeed, a number of experiments involving 
(p, n) reactions |I[ show that the additional channels of delta decay seem to dominate over 
the usual 7rN decay. 

Recently, we have proposed that the spreading width could be the source of errors in 
measuring the nucleon axial form factor from charged-current neutrino-nucleus scattering || . 
Namely, in the BNL experiment of neutrino-nucleus scattering || where neutrino beams with 
an average energy of 1.3 GeV are used, the data are usually fitted to the theoretical curve 
calculated in a Fermi gas model of the nucleus in order to extract the axial form factor. In 
this approach A-hole excitations are assumed to be completely excluded by experimentally 
rejecting events involving pions in the final state. However, when the A obtains a non-zero 
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spreading width, some events which generate A-hole excitations lead to pionless decays of 
the A which cannot be excluded by merely detecting the pions. To estimate these events 
quantitatively, it is important to know the branching ratio of the spreading width to the 7rN 
decay width before one can deduce the nucleon axial form factor accurately. 

Further, the knowledge of the branching ratio can be used to improve Monte Carlo 
simulations of neutrino detectors, e.g. in the Kamiokande experiment, where the flavor ratio 
of atmospheric neutrinos are measured |Q. It was suggested that pions produced from the 
charged- current neutrino-nucleus reactions cause uncertainties in identifying the flavor of 
the neutrinos created in the earth's atmosphere In these experiments the main source 
of pions is believed to come from delta decays. Therefore, a thorough knowledge of the 
branching ratio between pionic and non-pionic decays of deltas is crucial for estimating the 
pion events in neutrino-nucleus scattering of atmospheric neutrinos. 

In this work, we present a calculation of the delta decay width and its modification in 
nuclear matter within the framework of the Walecka model ||. The model is a relativistic 
meson-nucleon model which has been quite successful in describing nuclear properties, in- 
cluding response functions J/J. Crucial ingredients of this model are the presence of large 
vector and scalar self-energies which modify the nucleon properties in the medium enhancing 
relativistic effects. Our concern is to study the modification of the delta decay width driven 
by these relativistic effects. Nonrelativistic calculations of the delta decay in the medium 
can be found in Refs. ||U . 



II. FORMALISM 

In this section, we present the relativistic calculation of the A decay in the nuclear 



medium. First we review the calculation of the 7rN decay which has been studied in Ref. [ID 
Then we will present our calculation of the spreading width and its modification due to 
nuclear environment. 

In free space, a delta decays to 7rN. The decay width for this mode can be obtained from 
the imaginary part of the delta self-energy as shown in Fig. [I] (a). Formally, the decay width 
is given in terms of the self-energy E M1/ and the Rarita-Schwinger spinor of the delta, A M , 

r r -2Im(AV) (1) 

To calculate the self-energy, we introduce the 7rNA interaction Lagrangian, 

C«ka = — A M T N d^7r + h.c. . (2) 

Here N, 7r, and A are the nucleon, pion and delta field, respectively. The 2x4 isospin 
matrices T satisfy |L1 



T(T ] y = 5 ij - -tV . (3) 
3 

There is some ambiguity in the form of the Lagrangian of a spin-3/2 field, which affects 
the off-shell behavior of the particle. We will not discuss this complication in the following 
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because the delta is on-shell in our calculation. Using this Lagrangian and standard Feynman 
rules, one writes the self-energy in free space, 

= * J -0^%Q»D o (q)G f (k) , (4) 
where the pion [-D (g)] and nucleon [Gf(k)] propagators in free space are given as 

A>(?) = - — \— ~ > ( 5 ) 

<W = ■ < 6 > 

Putting these into Eq. ([[]), one can easily calculate the free decay width, 

r/= (M 2 %±^| kr m 



with 



■m v I 127rmA 



ml + m 2 — m 2 



k 2 = [ n n I -ml; £ k = Vk 2 + m 2 



2mA 



Here m, m^, and m^ represent the mass of nucleon, delta and pion respectively. By com- 
paring Eq. (^) with its experimental value of 115 MeV, one can fix the coupling constant 
/a = 2.15. 

In the nuclear medium viewed as a free Fermi gas of nucleons, the 7rN decay is suppressed 
by Pauli blocking, since the phase space available to the decaying nucleon is restricted by 
the occupied Fermi sea. Indeed, at nuclear saturation density the 7rN decay is completely 
blocked for a delta at rest in the medium. When the delta is moving, however, the decay 
sphere of the nucleon in momentum space becomes an ellipsoid, which only partially overlaps 
with the Fermi sphere ||10|| . Therefore, the delta obtains a non-zero energy- dependent width. 

The Pauli blocking can be incorporated by replacing Gf(k) with 

G {k) = G f (k) + 2ir5{k 2 - m 2 )Q{k )Q{k F - |k|)(# + m) . (8) 

Here kp is the Fermi momentum whose value at the nuclear saturation is 1.3 fm -1 . Note 
that in this form the nucleon propagator is given in the rest frame of the nuclear matter 
where the delta is moving with a certain relative velocity. From the following identity [TT 
for the A spinor with momentum p M = Pa) 



_ $a + rriA 



2m A ~ 3/2 

the decay width averaged over the delta spin becomes 

1 f JaV f d 4 k 



3 m A 3mA 3 

p^; 2 , (9) 



f.NdPAi) = ~ ^ / 7^e(go)% 2 - m 2 )6(fc° - E F )5(e - m 2 ) 

xTr[g^(^ + m)(^ A +m A )P^ 2 ] . (10) 
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Pa is the delta momentum relative to the nuclear matter, corresponding to a relative velocity 
v = J^-. Note that we have made use of Cutkosky rule [12[] in deriving this. The integration 
over the nucleon momentum k leads to 

2 
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i?± = max< E^, 
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A similar calculation has also been done in Ref . |TI| . The authors of Ref . [|1(J also discussed 
the results for different magnetic quantum numbers of the delta, which is qualitatively similar 
to the averaged result. 

The free Fermi gas calculation can be improved by incorporating mesonic mean fields. 
In a mean-field approximation to the Walecka model, the propagation of the delta and the 
nucleon is modified by the presence of constant scalar and vector mean fields which are 
determined from nuclear saturation properties. 

The 7rN decay width in this approximation can be obtained with the replacements 



m A - 
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(12) 



The nucleon self-energies (V^, Sn) are assumed to be the same as the delta self-energies (V^, 
5a) as is motivated from the universal coupling assumption. In symmetric nuclei, the space 
components of the vector self-energy are averaged to zero. Also its time component can be 
eliminated by appropriately shifting integration variables. However, the scalar self-energy 
Sn does contribute. Its value at a given nuclear density can be determined by solving the 
self-consistency equation 
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m 



(13) 



where the ratio of the scalar meson coupling to its mass, g 2 s /m 2 s , is determined from the 
equilibrium properties of nuclear matter 0. At nuclear saturation of ky — 1.3 fm _1 , Eq. ([13]) 
yields m* = 510 MeV. 

In addition to the 7rN decay, a delta in the medium can also decay into 2p-lh through 
the process, A + N — ► N + N. The lowest-order self-energy diagram for this process is shown 
in Fig. Q (b). This spreading width can be calculated in the same way as the 7rN decay 
by replacing the free pion propagator with the pion propagator D(q) including particle- hole 
polarization: 
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D(q) 



q 2 — ml — 11(g) + ie 



(14) 



We proceed by introducing the interaction Lagrangian for nucleons and pion using the 
pseudovector coupling for the pion, 



C 



7rNN — i " / lb i i " ^u,' 



(15) 



with /jy = 1.01 . Using the interaction Lagrangian Eq. ( |T5| ) we can easily determine the 
lowest order particle- hole polarization II ph (g): 



Re n ph (g) 
Im II ph (g) 



2 T7?, 2 f^F 




dE^ log 



(g 2 -2|k||q|) 2 -4g 2 £ 2 



(g 2 + 2|k||q|) 2 -4g 2 £ 2 



with 



E„ 



max 



Ei = Ep — E_ 
m, Ep - g , 



E_ = mm(E F , E max ) 




(16) 
(17) 

(18) 
(19) 



The A decay calculated with Eq. (14) contains two decay channels. First, the contribution 
from the usual irN decay which can be identified with the pole contribution of the pion 
propagator Eq. ([5]) with the shifted mass 



ml 



ml + n(g 2 



m 2 ' 



(20) 



If the lowest-order nucleon-hole polarization is used for 11(g), then, in some kinematics, 
Eq. (EO) is solved by negative m 2 indicating a tachyonic pole. This pole can be removed by 



introducing the Landau-Migdal parameter g' which is incorporated by writing [13] 



n(g) 



g 2 n ph 
g 2 + ^n ph 



(21) 



In our calculation, we take a conventional value of the Landau-Migdal parameter, g' = 0.7. 
One has to be a little bit careful with the pole contribution since the residue of the propagator 
is now changed due to the real part of 11(g). However, the modification of 7rN decay due to 
the real part of Eq. (^T|) is very small. 

The second contribution arises from the imaginary part of the polarization Eq. (|2~TD . This 
part can be identified with the particle-hole decay channel of the virtual pion generating the 
spreading width of the A. Note that the imaginary part of the particle-hole polarization 
is only non-zero at space-like momenta of the pion. Therefore, the spreading width can 
be separated from the 7rN decay in the medium due to their different kinematic region. 
Replacing D (q) by the expression Eq. ( |T4"[ ) in Eq. (^) and taking the contribution arising 
from the imaginary part of the polarization, we obtain the spreading width in the nuclear 
rest frame, 
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f.(|PA|) = 

[g 2 -ml- ReII(g)P + [ImH(g)P ^ + ™ m ^ I* m * ~ ' «) \ > ( 22 ) 

which can be solved numerically. The lowest order spreading width shown in Fig. (|I|) (b), 
which is usually referred to as the spreading width, can be obtained from Eq. (p^ ) by taking 
the lowest-order contribution from the polarization, only, and by neglecting the Landau- 
Migdal correction. Note that, unlike in the case of the free 7rN decay, the pion momentum, 
q 2 , is not restricted to its on-shell value in the course of the integration. Therefore, one 
should consider the momentum dependence of the 7rNA and 7rNN vertices. In this regard, 
we employ the monopole form for the form factor, 

F ^ = • < 23 > 

in our discussion. This is to take into account the finite size of the nucleons as well as the 
short-range effects of NN interactions. However, the value of the cut-off parameter A is 
in doubt. In principle, this should be determined by summing the vertex diagrams of the 



underlying field theory. Here we use the Bonn potential values ||14| ; for 7rNN, A = 1.8 GeV 
and for 7rNA, A = 0.85 GeV. Finally, in order to implement the effects due to the nuclear 
self-energies, the replacements in Eq. (|12|) are done for the calculation of the spreading 
width. 

In addition to the particle-hole polarization, A-hole excitations can contribute. This 
means that a delta rescatters in the medium by A + N — ► A + N, which can be obtained 
by taking the imaginary part of A-h polarization. The difference to the particle-hole polar- 
ization is that the imaginary part has finite values even for positive q 2 . Depending on the 
four-momentum the pion pole is shifted into the complex plane as the pion obtains a width 
from the decay channel to a A-h excitation in the medium. The rescattered A subsequently 
decays leading to more complex decay channels, for example 2N - lh - In, etc. The possible 
influence of this contribution will be studied in a forthcoming publication ||15|| . 



III. RESULT 

We now present the result for the various delta widths. Our concern is to understand 
the momentum dependence of the decay width as well as the role of the nuclear mean fields. 

In the following, we will show our results calculated at half, normal and twice the nuclear 
saturation density with the corresponding Fermi momenta kp = 203 MeV, 256 MeV, and 
323 MeV, respectively. The self-consistency equation of Eq. (|T3|) yields an effective nucleon 
mass m* = 510 MeV at the saturation density, m* = 715 MeV at half saturation density, and 
m* = 261 MeV at twice nuclear density. The reduction of the nucleon mass in the medium 
arises due to the strong scalar self-energy. In principle, the self-energies are functions of 
the momentum of the propagating nucleon [l(J. Specifically, the absolute values of the 
self-energies decrease as the momentum increases. This momentum dependence could be 
interesting especially in the high momentum region. However, in most of the momentum 
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range of our concern, the momentum dependence of the self-energies is not expected to be 
important. The result with half the saturation density is applicable to (p, n) reactions where 
the interactions most likely occur on the nuclear surface. On the other hand, the result 
at saturation density could be used in lepton-nucleus scattering as the incoming lepton 
can probe the region deeply inside the nucleus. The result with higher density could be 
interesting in the heavy-ion collision even though the neglect of contributions from the lower 
continuum of the nucleons in our calculation is questionable. 

Figures ^| show the results for the 7rN decay width as a function of the delta momentum 
| Pa |- Fig. |2| (a) shows the result of free Fermi gas calculation while Fig. |2| (b) includes 
the nuclear mean fields. At saturation density or at twice the saturation density, the decay 
width increases as the delta momentum increases. This is because the overlap between the 
momentum ellipsoid of the decay nucleon and the Fermi sphere is getting smaller and the 
Pauli blocking becomes less important. Indeed, as is shown by the dashed curve, the ellipsoid 
is completely out of the Fermi sphere for |pa| > 650 MeV, and the 7rN decay width is the 
same as its value in free space. A similar trend can be observed for the mean field calculation 
in Fig. |2| (b). At half saturation density (the solid curve in both figures), however, Pauli 
blocking is less significant as could be expected. Our results indicate that the mean field 
combined with the Pauli blocking reduces the 7rN decay width substantially from its free 



value of 115 MeV which agrees with Wehrberger and Wittman [|10 . 

Unlike the 7rN decay, the spreading width, which is represented by the mechanism A + 
N — > N + N, is plotted in Figs. [| Figure |3| (a) shows our calculated spreading width without 
the mean fields and Figure |3| (b) includes the mean fields. First note that the momentum 
dependence of the spreading width is not as significant as the 7rN decay. This means, when 
added to the 7rN decay width, the momentum dependence of the in-medium decay width is 
mainly driven by the Pauli blocking. 

Conventionally the spreading width is scaled as the nuclear density ]17| . This is reflected 
in the Fermi gas results in Fig. |3| (a) where the dashed curve is about a factor of 2 smaller 
than the dot-dashed curve and the solid curve is about a half of the dashed curve. But this 
kind of scaling is not applicable in the mean field case. Therefore, scaling the spreading 
width with density has to be applied with caution. Our mean field results at the saturation 
density which vary from 75 MeV to 57 MeV qualitatively agrees with the conventionally 



adopted value of 70 MeV in the nonrelativistic calculations fT7j , [T8|j . Also at this density, 



the branching ratio of pionic decay to the spreading width, T^/T s , varies from to 1.33 
indicating that a significant portion of deltas leads to the pionless decays. However the role 
of the spreading width compared to the 7rN decay becomes smaller at half saturation density. 

On the other hand, the mean field result for the spreading width at the double saturation 
density is rather abrupt. Its magnitude is only 6 ~ 12 % of the Fermi gas result. This kind of 
huge reduction due to the mean fields might come from neglecting the vacuum contribution 
to the particle-hole loop in Hartree approximation. At this density in our simple calculation, 
the effective mass of nucleon is 261 MeV, much smaller than its free value. Therefore, anti- 
nucleons can be easily excited which could lead to the appreciable contribution from the 
vacuum. More work on the vacuum contribution should be done ||15|| . However, as indicated 
in Fig. |3] (a) and (b), the spreading width has a very strong dependence on the mean field 
given our approximations. 

A-hole contributions in neutrino-nucleus scatterings have been discussed in Ref. [0,§J. 
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There it was stressed that the pionic versus non-pionic decay could be an important source 
of uncertainty in determining the nucleon axial form factor or estimating the pion events in 
atmospheric neutrinos. For example, in Ref. 0, we have found that the theoretical error 
of about 0.1 GeV in the extracted axial mass M4 is expected from the events due to delta- 
hole excitations. This conclusion is drawn by assuming that the branch ratio between the 
spreading width and the 7rN decay is roughly one. However, our calculation in this work 
show that the branching ratio has a clear momentum dependence which should be properly 
taken into account. In this regard, our results, especially the results at saturation density, 
could provide important constraints. 

In summary, we have calculated the delta spreading width in the framework of a rela- 
tivistic nucleon-meson model. Our concern is to present the branching ratio of the spreading 
width to 7rN decay so that it can be used for the analysis of measurements of the nucleon 
axial form factor || or for the atmospheric neutrino experiments [|J. 7rN decay is substan- 
tially reduced by Pauli blocking and mean fields. Also, mean fields play an important role 
in the spreading width. We have found that the spreading width is large even though its 
relative strength to the 7rN decay depends on the nuclear density and the relative motion of 
the delta with respect to the nuclear medium. 
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FIGURES 



FIG. 1. Diagrams for A self-energy. Figure (a) is the self-energy diagram representing the 7rN 
decay and figure (b) is the lowest diagram for the spreading width. 

FIG. 2. Pionic decay of A with respect to the delta momentum, (a) is for the 7rN decay without 
mean fields and (b) includes the mean fields. The solid lines show the results for the half saturation 
density, the dashed lines for the saturation density and the dot-dashed lines for twice the saturation 
density. 

FIG. 3. Spreading width of A. (a) shows the result without mean fields and (b) includes the 
mean fields. The solid lines show the results for the half saturation density, the dashed lines for 
the saturation density and the dot-dashed lines for twice the saturation density. 



10 




Figure 1 



11 




Figure 2 
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